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Abstract 

Water is an associated liquid in which the main intermolecular interaction is 
the hydrogen bond (HB) which is hmited to four per atom, independently of 
the number of neighbours. We have considered a hydrogen bond net super- 
posed on Bernal's geometric model for liquids, which allows for different local 
environments for the liquid particles. In this study, a mean-field treatment 
of the two-dimensinal version of the model is discussed. Under pressure the 
model exhibits three phases of different densities and a coexistence line ending 
in a critical point between low and high density phases. Entropy of the HB 
network plays an essential role in defining the slope of the coexistence line. 
The model behaviour might be of interest in describing supercooled water and 
liquid-liquid transitions of other substances. 
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I. INTRODUCTION 



Liquid polymorphism has been proposed for different systems, both from experimental 
evidence, as well as from theoretical arguments |l|-10|. The relevant systems must have open 
molecular coordination environments at low pressures, and may either be liquids associated 
through a hydrogen-bond net, water being the main example or substances for which the 
different possibilities for hybridyzed eletronic orbitals imply structures of different densities, 
as in the case of carbon 

The anomalous behaviour of water has intrigued physicists and chemists for a long time 
TTHT4|. It is presently accepted that the special properties are related to hydrogen bonds. 



Because the energy involved in a hydrogen bond is an order of magnitude larger than the 
typical van der Waals energy fT^JTBH , the fusion and boiling temperatures of water are much 
higher than those of homologous substances [jl4|. In fact, liquid water exhibits an extense 
random network of hydrogen bonds which continuously reformulates, but presents local 
tetrahedral symmetry, responsible for the large volume at low temperatures • 



In the last years, interest in supercooled water arose P,p!7|-p!9[| as a means to explore the 
behaviour of isothermal compressibility and specific heat which, differently from most sub- 
stances, rise as one lowers temperature towards the fusion point and beyond p^j20[| . Search 
for an explanation of these two properties led to three different hypotheses: a retracing 
liquid-gas spinodal a second critical point in the supercooled region P JT7|j2l|| , and a 
"singularity-free" scenario p2|- p^] , in opposition to the first two. Sastry and collaborators 
have shown that a singularity is not the only explanation for the compressibility or 
specific heat low temperature behaviour, but molecular dynamics simulations of ST2 wa- 
ter p5| indicate that a divergence of compressibility could occur around 2000atm, 230K . 



Among the several simplified models presented in the literature, some do p6|-p8|, other do 



not p2|-p^ present the second critical point, in some cases depending on model parameters 
used. From the experimental point of view, the question is very difficult to settle, because 
in the supercooled region one encounters a nucleation line (— 38°C, at 1 atm) beyond which 
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the system spontaneously crystallizes 0. 

The second critical point of water would be the end of a coexistence line between two 
metastable liquids with different densities (HDL and LDL for high and low density liquid) 0. 
These two liquids would be related to two amorphous solids which exhibit coexistence and 
are found experimentally beyond the nucleation line . The existence of the amorphous 



phases lends support to the idea of a liquid-liquid transition. 

The idea of liquid polymorphism has been applied not only to water, but also to certain 
molecular liquids, such as I, Se and S, in order to explain an abrupt increase in conductiv- 
ity, which would be related to a first order phase transition in the supercooled region [|l|. 
Experimental evidence of a transition between two distinct thermodynamically stable liquid 
phases in carbon was found [Q] some years after a suggestion of its existence appeared in the 
literature [Q. This transition has also been reported from molecular dynamics simulations 
and a scaling formalism ||^. Experimental results also report a liquid-liquid transition in 
phosphorus in the stable phase 0]. 

A few decades ago Bernal proposed a model for liquids which consisted of a close-packed 
assembly of different regular or quasi-regular polyhedra . The typical hquid configuration 
was expected to be statistically homogeneous and to possess no long-range order, due to the 
presence of different kinds of polyhedra. Molecules would be located only at the vertices 
and the edges would represent an average nearest neighbour intermolecular distance. The 
absence of holes, one of the model hypotheses, reduces the number of allowed polyhedra. 
Collins studied a mean- field solution for a version of Bernal's model in two dimensions ||31|| , 



in which case the polygons to be considered are squares and triangles. Collins demonstrated 
that the model could exhibit a phase transition, at fixed volume, for a specially designed 
set of interaction constants with no physical motivation. The author was looking for a 
melting transition, but recognized that the phase change was more analogous to a change of 
association number from one liquid to another. Interestingly, a recent study ||32| of a dense 



system of repulsive Lennard- Jones particles showed that the spatial configurations can be 
interpreted in terms of random square-triangle tilings, as in Bernal's proposal. A transition 



between two densities was also found. 

We have considered the spatial distribution of particles proposed by Collins and intro- 
duced interactions between the molecules which would mimic the HB network. Thus the 
distribution of energies is a result of the disponibilities of molecules to engage in HBs and 
favours low-density arrangements at low temperatures and pressures. Also, besides the 
translational entropy for simple liquids, an entropy term related to the distribution of HBs 
must be considered. The model does not allow for vacancies, so a gaseous phase is dis- 
regarded. A mean-field treatment of the system behaviour as a function of pressure and 
temperature is undertaken. We have also considered a variant of the model, which could 
describe non-associated liquids with preference for low density at low pressures (that would 
be the case of carbon). In the latter case, the HB entropy contribution is absent. 

Polymorphism and a critical point are found for both models. However, the slope of the 
coexistence lines depends on competing entropy contributions. The model for water and 
the corresponding phase diagram is presented in section II. Liquid-liquid transitions and 
the phase diagram for non-associated liquids are discussed in section HI. Discussion and 
summary are presented in section IV. 

II. THE MODEL FOR WATER 

A. The geometric description 

In two dimensions, the Bernal liquid consists of a system of adjacent and randomly 
distributed squares and triangles of equal sides. The particles are localized at the vertices 
and there are three possibilities for the number of nearest neighbours r, as showed in Fig. 
|l]-a. A 4 and a 6- molecule can never be nearest neighbours by a geometric constraint. 
We call rij. the number fraction of molecules with coordination r (r-molecule) and v^ the 
corresponding specific volume, 
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(1) 



for fixed intermolecular lenth h. 

Writing v = V/N for the volume per molecule, number and volume conservation are 
written as: 

E ^r = l (2) 

r=4,5,6 

UrVr = V. (3) 

r=4,5,6 

Interactions must still be specified. In simple liquids, van der Waals interactions depend 
mainly on interatomic spacings. Hydrogen bonding depends on disponibility of neighbours 
to accept or donate hydrogens, thus in liquid water loss of translational order and increased 
density results in the frustration of some of the HBs for molecules with over four near- 



est neighbours In order to represent this property we introduce, alongside with the 

distribution of triangles and squares, a distribution of bonds. We have considered four as 
the maximum number of bonds per molecule [^], thus the HB between a 4 — 5 pair, for 



instance, may be absent, depending on the distribution of bonds among the neighbours of 
the 5-molecule (see Fig. l-b). Note that the energy of the pair depends on the distribution 
of HBs amongst their neighbours. Because the hydrogen bond is an order of magnitude 
larger than van der Waals interaction [|12|, in this study we have considered only the HB 
interaction. The directionality of the hydrogen bond was also ignored. 



B. A mean- field approach for the HB network 

Our approximation for the distribution of HBs is the following: the four possible bonds of 
an r-molecule are distributed randomly over the r possible lines and the average energy of an 
r — s pair is calculated as resulting from the two independent distributions (for example: for 
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5-molecules there are five possible distributions of four bonds on five lines, so the probability 
of an HB for a 5 — 5 pair is ^). In this manner, if the HB energy is — e, then the average 
r — s energy per bond is given by: = which makes $44 (the energy of a pair of 

4-coordinated particles) the minimum energy. Under this assumption, the energy of the 
model may be written as: E = J2rs ^rs^rs where Nrs is the number of rs pairs. 

Standard mean- field conditions may then be implemented by assuming N^s = rNrPs{r) + 
sNsPr{s), where Nr = rirN is the total number of r- molecules and Ps{r) is the probability 
that an r-molecule has an s-molecule as neighbour. Notice that the geometric restriction 
must be taken into account, so assuming random distribution of molecules we have: 

= A 

where Ntot{r) is the maximum number of neigbours for an r-molecule {Ntot{4:) = N4 + N^, 
Ntot{5) =N^ + N, + Nq and Ntot{6) = + Nq). 
The average energy per molecule is then: 

e{{ni}) = ^J2^^rJ2Psi^)'^rs, (5) 
^ r s 

where | stands for double counting. 

C. Entropy and Gibbs free-energy 

The partion function of the system in the Gibbs (constant pressure) ensemble is written 

as: 

Z(r,P) = ^fi({n.})e-^, (6) 

{rii} 

where Y^' represents the summation over rij constrained by particle number conservation 
(Eq. 0), and Q is the number of states with volume V{{ni}) = Nv{{ni}) and energy 
E{{ni}) = Ne{{ni}) (Eqs. |^ and Two factors contribute to entropy, one related with 
the spatial distribuition of particles, ^partiini}) = ^vlL^Ovel' other related with 



6 



the distribuition of HBs, flnBiini}). The degeneracy of HBs is taken into account in the 
simplest way, by considering independent molecules: 

nHB = l[D:^^ (7) 

r 

where Dr = ^r^^^L^y is the number of possible arrangements of four HBs over the vertices of 
an r-molecule. Thus Z (Eq. |^) is rewritten as: 

. ' . JV7({n, }) 

Z(T,P) = Ee-^^, (8) 

where A^7 = E — TS + PV and 5* = klnfltot {^tot = ^part^Hs)- In the thermody- 
namic limit the Gibbs free energy per molecule is determined by the equation: g{T, P) = 
"^{T, P; {njjmin)) where {njjmin is the set that minimizes 7. In order to find it we minimized 
7 with the constraint of Eq. ^ obtaining two transcendental equations, which are solved 
numerically. Henceforth the temperature and pressure will be written in reduced units t and 
p {KBT/e and Pb'^/e, respectivelly). 

D. The phase-diagram 

We find the equilibrium values of function of temperature and pressure using a 

combination of numerical methods. A rough determination of the free energy minima in an 
extensive search was followed by the Newton method, for more precision. 

The phase diagram of the model for liquid water on the pressure-temperature plane is 
shown in Fig. |^. At lower temperatures, the model presents two coexistence lines between 
three phases with different densities, with a triple point at an intermediate temperature. 
At higher temperatures, the coexistence between a high density (HD) and a low density 
(LD) phase ends at a critical point. The critical temperature is obtained both from the 
study of the free-energy surface as well from the evolution of the order parameter defined as 
A = Phd — Pld where Phd,ld are the corresponding densities. 

The triple and critical points (in reduced units) are given respectivelly by: ttp = 0.0799, 
Ptp = 4.0780, tc = 0.0965 and pc = 3.9241. Curiously, if one assumes 20kJ/mol [|lll for the 
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HB energy (e), a critical temperature of 232K is obtained (against 230K obtained from MD 
simulations [0). As for the critical pressure, dimensionality would enter the calculation, 
rendering it meaningless. 

Fig. ^a illustrates discontinuity of volume on the coexistence lines. The low density 
phase has more four-coordinated molecules at low temperatures, but on the coexistence line 
between LD and HD six-coordinated particles predominate in both phases, as can be seen 
in Fig. ^b. The three phases differ basically in density. 

The negative sign of the slope between the HD and LD phases in Fig. ^ indicates, in 
accordance with Le Chatelier's principle, that entropy increases while volume decreases, on 
transition into the high temperature phase (HD), as can be seen in Figs, ^a and |^c. The 
higher density phase has higher entropy than the LD phase. Although it is true that the 
completely ordered 6-phase has higher entropy than the completely ordered 4-phase, because 
of entropy of bonds (see Eq. |^, at finite temperatures competition with translational entropy 
may invert the balance in some regions of the phase diagram, as can be seen in Fig. ^ at 
very small temperatures. But the overall picture is that the higher density phases also have 
higher entropy. 

Compressibility and specific heat for the model can be obtained numerically. These quan- 
tities present maxima whose magnitude increase with decreasing temperature, as expected 
on approach of the critical point. The line of specific heat and compressibility maxima is 
shown in Fig. |. The distance between the two lines grows as pressure is lowered. The 
specific heat maxima are localized at lower temperatures with respect to the compressibil- 
ity ones, a feature also found for the model discussed by Sastry and collaborators, which 



presents no liquid- liquid transition [22 



III. LIQUID-LIQUID TRANSITIONS 

In case one does not consider the HB net entropy (Eq. 0), a different phase diagram is 
obtained. As seen in Fig. R, the intermediate density phase disappears and an inversion 
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occurs with respect to entropy: the higher entropy phase, in this case, is the low density 
one, as can be guessed from the slope of the coexistence curve. These features might be 
related to other liquid-liquid transition such as that of carbon. 

The different structures of solid carbon, diamond and graphite, are a result of hybridiza- 
tion of the carbon atoms. The diamond structure, stable at higher pressures, is associated to 
sp^ tetrahedral hybridization, while graphite is constituted of planes of sp"^ hybrids. Graphite 
is less dense, but of higher entropy. Experimental and theoretical results point to the 

possibility of a liquid-liquid phase transition, between a high density liquid phase (dominated 
by sp^ hybridized atoms) and a low density liquid phase (dominated by sp hybridized atoms 
in the case of MD simulations 0] and by sp"^ hybridized atoms in a scaling formalism model 



gl). As in the sohd state [g4|, the slope of the liquid-liquid coexistence hne in the pressure 
versus temperature diagram is positive and, according to Le Chatelier's principle, the 
lower density phase must present the largest entropy. In the absence of the HB net entropy, 
the model we propose could be thought of as a mean-field treatment of Brenner's potential, 
used to describe carbon, and for which bond energies depend on the local environment in 
such a way as to produce the correct geometries and energies of the known carbon structures 
PJ53[|. Is is not clear whether it would be necessary to consider additional entropy terms. 



because of the mean-field treatment. As a curiosity, in our model, from tc = 0.062 (Fig. 
if we consider carbon bond energy of the order of 600k J /mol we obtain ~ 5000A' 
(against 9000/^ from MD simulations 0). 

The two cases studied lead to a simple description of potential candidates for liquid 
polymorphism. We discuss this in terms of ananalysis of minimization of the Gibbs free- 
energy. At small temperatures, at which energy wins out against entropy, pressure may 
produce a transition betweeen a high and a low density 2/ the low energy phase is the less 
dense (AG ~sma//T AU — P\AV\ < at some pressure, where AX is defined as Xhd — ^ld)- 

Thus liquid-liquid transitions may arise if low lying energy states are low density states 
(as is the case of water or carbon) and for any model with energy monotonically increasing 
with density. Of course, in order for the system to present also a gaseous phase, mean field 



or van der Waals energy should increase towards very low densities, and thus present an 
extremum, as pointed out by Tejero and Bans ||T0| and also implied in Poole and collabora- 
tors' proposal ||28|. However, we understand the van der Waals attractive potential would 
have to present a minimum (as can be inferred from I^Hj): as density increases from the gas 



region, the potential decreases, as in usual van der Waals theory, but a minimum is present 
at some intermediate density, above which the potential slightly increases, destabilizing a 
very dense liquid phase, as in this model. Once the condition for the liquid-liquid transition 
is fulfilled, two possibilities arise for the slope of the coexistence line, and these, according 
to Le Chatelier's principle, will depend entirely on entropy: if the high density phase is the 
higher entropy phase, as in water, we encounter a negative slope, whereas if it is a smaller 
entropy phase, as in carbon, we encounter a positive slope. 

IV. DISCUSSION AND SUMMARY 

In liquid water there are two sources of disorder, not really independent: positional 
disorder of the molecules and disorder of hydrogen bonds. In the case of the latter, molecules 
may loose HBs either (i) because there are no nearby molecules free to engage in an HE (as 
in the case of a disfavourable increased local coordination - see Fig. |l]-b) or (ii) because of 
rotations which remove the necessary alignment between neighbouring molecules. In some of 
the models proposed in the literature only the second feature has been considered [^p8 



m 



terms of ice-like models, while in others p2| , |23| , |26| molecular rotations leading to weak bonds 
are described through Potts variables on cubic lattices, and smaller volume attributed to the 
disordered bonds in order to account for the first feature (increased non-alignment in dense 
environments). Interestingly, both kinds of model (except p6[) present only a liquid-gas 
transition and no low-temperature transition, in support of the "singularity-free" scenario. 
A low-temperature transition in the ice model would imply a discontinuity in the number of 
HBs at the transition, a requirement which would also have to be met in the Potts models 
(except [^), for which volume is measured directly in terms of bonds. 
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In the case of the model discussed in this paper, only the first feature, the distortion of 
the HB net, leading to rupture of HBs, not due to non-alignment, but to disponibility of 
receptor or donor molecules, is taken into account. In this case, liquid-liquid transitions are 
present, independently of model parameters, but for specified geometric constraints for local 
structure. 

In relation to the dispute on the existence of two supercooled waters , the answer could 
depend on which of the two mechanisms above, restricted bonding due to varying local 
coordination, or weak bonding due to thermally induced non-alignment, is most relevant. 

Our study of the Bernal model for liquids on the pressure-temperature plane has shown 
the possibility of liquid-liquid coexistence if energy favours low density configurations. Liquid 
polymorphism both above and below the fusion line have been suggested in the literature 
for associating, as well as for non-associating liquids. Our results point to the role of bond 
network entropy. The latter becomes patent if one compares Figs. ^ and |^. The inclination 
of the coexistence lines present opposite signs in the critical point region, for the two systems. 

The compressibility and specific heat maxima beyond the critical temperature could 
relate to the so called anomalous properties of water. However, as for other models discussed 
in the literature, it is senseless to discuss whether the polymorphisms belongs to supercooled 
or stable liquid phase, since it is unable to describe the solid phase. Also, the model does 
not allow for "holes" and is therefore uncapable of describing the liquid-gas transition, but 
this is not really a shortcoming, because interest lies really in the dense phases. 
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FIG. 1. (a) Three possibilities for the local molecular environment: 4, 5 and 6-molecules (black dots), 
(b) The existence of an HB between a 4 and a 5 molecule depends on the distribution of the four molecular 
bonds (full line) over the respective "neighbouring" lines (dashed lines). In the left figure no HB occurs 
between the two central molecules because the 5-molecule has engaged its four bonds with the other 4 
neighbours. 
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FIG. 2. Coexistence lines and pressure-temperature phase diagram in reduced variables. Three phases 
of low, intermediate and high density are identified (LD, ID and HD). The LD phase has n4 > hq > 
for low temperatures, uq > n4 > for t > 0.06 and ne > ns > m above t = 0.09. ID and HD have, 
respectively uq > n4 > and ne > > n4, independently of temperature. The triple and the critical 
points are indicated. The inset shows the region of interest for supercooled water. 
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FIG. 3. (a) Volume versus temperature on the coexistence lines. Coexistence regions are indicated by 
gray lines. LD, ID and HD as in Fig. 2. (b) n4 (full lines), (dashed lines) and uq (long dashed lines) 
versus temperature for the LD (circles) and HD (squares) phases on the coexistence line above the triple 
point. For the region of temperatures shown six-coordinated particles predominate in both phases, (c) 
Entropy of the HD (squares) and LD (circles) phases versus temperature on the coexistence line. 
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FIG. 4. Lines of specific heat (squares) and compresibility (circles) maxima and coexistence line (tri- 
angles) near the critical point. 
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FIG. 5. Coexistence line on the pressure-temperature phase diagram in the absence of HB net entropy. 
The high density phase (HD) has uq > Ui and the low density phase (LD) presents Ui > ne for low 
temperatures and ng > n4 for t > 0.06. The critical point is given by tc = 0.0619 and Pc = 5.2843. The 
inset shows reentrant behaviour at very low temperatures. 
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